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SHARP INTERFACE LIMIT OF AN ENERGY MODELLING 
NANOPARTICLE-POLYMER BLENDS 

STAN ALAMA, LIA BRONSARD, AND IHSAN TOPALOGLU 


Abstract. We identify the P-limit of a nanoparticle-polymer model as the number of particles 
goes to infinity and as the size of the particles and the phase transition thickness of the 
polymer phases approach zero. The limiting energy consists of two terms: the perimeter of the 
interface separating the phases and a penalization term related to the density distribution of 
the infinitely many small nanoparticles. We prove that local minimizers of the limiting energy 
admit regular phase boundaries and derive necessary conditions of local minimality via the 
first variation. Finally we discuss possible critical and minimizing patterns in two dimensions 
and how these patterns vary from global minimizers of the purely local isoperimetric problem. 


1. Introduction 


In many applications, engineering a self-regulating, stable structure with predetermined 
physical properties is highly desirable. Here we consider the case of block copolymers, for 
which a composite is created by adding solid “filler” nanoparticles in a blend of macromolecules 
to create high-performance polymers that are used in, for example, solid-state rechargeable 
batteries, photonic band gap devices, etc. (cf. [3j and references therein). Depending on the 
desired physical properties, block copolymers can be used to direct assembly of nanoparticles 
or, vice-versa, nanoparticles can be placed in the polymeric matrix to alter the morphology 
of block copolymer microdomains in both the strong and intermediate segregation limits (see 
e.g. [33 37 38l|55]). Nanoparticles are also used in altering morphologies of immiscible mixtures 


such as oil/water mixtures or fluid-bicontinuous gels with a microreaction medium and in 


controlling domains of coarsening (cf. [16,21,38,53 


In this paper we consider an Ohta-Kawasaki-type model for a nanoparticle-block copoly¬ 
mer composite and study its limit as the interfacial length scale becomes small, within the 
framework of T-convergence. The limiting sharp-interface model combines three effects via 
interfacial, repulsive nonlocal and bulk terms. As a first step, in this paper we concentrate 
on two of these and provide a more detailed analysis of a sharp-interface energy consisting 
only of the interfacial and bulk terms. Clearly, retaining the repulsive interaction term would 
enrich the energy landscape of the problem and admit a wider variety of potential patterns 
for minimizers; however, by excluding nonlocal effects we may focus on the confining effect of 
the nanoparticles which already dramatically alters the morphology of the energy minimizers 
for the isoperimetric problem. We illustrate this with a specific example in two dimensions, in 
which the presence of nanoparticles influences the minimizing configuration to switch from a 
lamellar to a circular interface. 
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1.1. The Limiting Problem. We assume that the nanoparticles have a given, fixed distri¬ 
bution, described by an absolutely continuous probability measure p £ & ac (T n ). The sharp 
interface limit is to minimize the energy 

(1.1) E ^ a (u)\=\f \Vu\+a [ (u(x) - l) 2 dp(x). 

Z Jjn Jfn 

over functions u £ BV(T n ; {±1}) satisfying a mass constraint 


( 1 . 2 ) 



u(x) dx = m 


where m £ (—1,1) and a > 0 are constants and T n is the n-dimensional flat torus. Here 
f Jn |Vit| denotes the total variation of the function u £ BV (T n ; {±1}) and is defined as 


f |Vii| := sup | f udivtpdx: tp £ C'o(T n ;M n ), |^>(x)| 

Jjn l J T n 


^ 1 


Hence, the first term in the energy calculates the perimeter of the interface between the phases 
1 and —1 whereas the second term penalizes the phase u = — 1 according to the given prob¬ 
ability measure /i. The strength of this penalization is controlled by the parameter a. The 
energy arises as the singular limit of a sequence of energies that appear in models of self- 
assembly of nanoparticle-polymer blends and are given by the standard Cahn-Hilliard energy 
with an inhomogeneous term. This limiting energy can be considered as the extension of the 
classical isoperimetric problem to an inhomogeneous medium. It also gives a simple model of 
understanding how penalization of one phase via a probability density affects the geometry of 
the phase boundary. Indeed, depending on the choice of the measure p (and the mass m,) 
the penalization term here can act as an attraction or repulsion between associated phases. 
Minimization of energies similar to E Mi(T where the competition between the interfacial and 
bulk energetic terms drives a pattern formation has attracted much interest in the past. In 
particular, energies with a penalization term in this spirit appear in modelling image process¬ 
ing problems such as image segmentation, inpainting and denoising (cf. (9j[l8j[3TJ[32]). This 
energy also has much in common with the problem of minimizing perimeter in the presence of 
an obstacle constraint (cf. Bill [22j[34]). 


1.2. The Diffuse Interface Model. To model the nanoparticle-block copolymer configura¬ 
tions with N nanoparticles where each particle is of the form of an n-dimensional ball of radius 
r > 0 and center x\ £ T n for i = 1,... ,1V, in 1 19 20 the authors consider a free energy which 


extends the Ohta-Kawasaki model 45 with an addition of a penalization term. In a dynami¬ 
cal model, one would expect the nanoparticles to be mobile as they interact with the diblock 
copolymers. Here we assume that the nanoparticle dynamics is at a significantly slower time 
scale; hence, we fix their location, and treat them as confining or pinning elements of polymer 
chains. Neglecting the mobility of nanoparticles is reasonble since we would like to observe 
the change in the polymer morphology in terms of the minimization of the free energy with 
respect to the phase parameter. To be precise, in its most general form we consider the free 
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energy 


(1.3) 


-e,^f,7],m,r,u p ,N 


(u;x):=^/ \Vu\ 2 dx + f {u 2 - l) 2 dx 
o J T n J T n 



G(x, y)(u(x) — m)(u(y) — m) dxdy 


V VV(| x - Xi\){u- Up) 2 dx. 

J'Tn 

J i=l 


Here T n denotes the n-dimensional flat torus, x = (aq, .. ., xjv) £ T n x • • • x T" denotes the 
vector consisting of the centers of nanoparticles and u E H l { T”) is the phase parameter. 
Moreover, as above m = f Tn u(x)dx, and e>0, 7 > 0 , 77 > 0, r > 0 and u p E [— 1 , 1 ] 
are constants. Clearly m E (—1,1) is related to the volume fraction of the diblock copolymers 
describing the distribution of different types polymers with m = 0 corresponding to equal mass 
distribution between two types of polymers. As it is standard in Cahn-Hilliard-type energies, 
e > 0 describes the thickness of the transition layer between u = 1 and u = — 1. (The extra 
factor of j will simplify the form of the eventual Gamma limit, and is inconsequential.) The 
parameter 7 > 0 is related to the strength of the chemical bond between the polymers subchains 
in the copolymer macromolecule and controls the long-range interaction between the phases 
via the Green’s function of the flat n-torus denoted by G(x,y). The parameter u p E [— 1 , 1 ] 
determines the “wetting” of nanoparticles and their preference towards the polymer phases. 
For example, u p = 1 means that the particles prefer to stick to those subchains of a diblock 
copolymer macromolecule that are given by the phase u = 1. Finally V denotes a rapidly 


decreasing repulsive potential. In 19,20 , the authors take V(|x|) = exp(—|x|/ro) with ro <C 1 
so that the repulsion is short-ranged; however, we restrict V to be a smooth and radial function 
of support B( 0, r ) so that the interaction between particles is zero. Hence, the last term of the 
energy simplifies to 


N 

E 

i= 1 ' 


B(xi,r) 


V (|x — Xi\)(u(x) — u p ) 2 dx. 


Remark 1.1. (The potential V) The reason for the inclusion of a potential V as a weight in 
the penalization term is related to the dynamics of the model. Indeed, in their model the 
authors consider the evolution of a system of nanoparticle-block copolymer blend as a gradient 


flow of the free energy given by (1.3) along with a system of ordinary differential equations for 


evolution of the centers of nanoparticles. There a nonzero potential V enables the nanoparticles 
to move around in the domain T n whereas its short-range allows the authors to neglect the 
interaction between nanoparticles. 

Even in the absence of nanoparticles (i.e., when rj = 0) the microphase separation of diblock 
copolymers yields a rather rich and complex picture. There is an extensive literature on the 
mathematical analysis of phase separation of block copolymers via the Ohta-Kawasaki model 
and its sharp interface limit leading to a nonlocal isoperimetric problem. From mathematical 
derivation of the model [8,12 to analysis on curved manifolds 15,56 , the energy landscape of 
(1.3) with i] = 0 and its T-limit as e —> 0 whether posed on the flat torus (i.e. with periodic 


boundary conditions), on a general domain with homogeneous Neumann data or on the whole 
Euclidean space has been rigorously investigated in various parameter regimes of m and 7 
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(cf. [lj[2j|6 10 11,25 26 35,36 39 42-44,46-48 51 ). However, to our knowledge, mathematical 
analysis of nanoparticle-block copolymer blends via the energy (1.3) or its sharp interface 


version ( 1 . 1 ) has not been carried out. 


1.3. Choices of parameters. We concentrate only on the local interactions between the 
phases, i.e., we choose 7 = 0. This simplification does not affect the passage to the T-limit in 
the full energy functional, as the nonlocal interaction term is a continuous perturbation of the 
other terms, but (as explained above) in this paper we restrict our attention to the competition 
between the interfacial and bulk confinement terms only. In the absence of nanoparticles 
(7 = 0) the periodic phase separation and the relation between the periodic Cahn-Hilliard 
energy and the periodic isoperimetric problem have been investigated in 0- The addition of 
nanoparticles into the model, of course, poses new challenges. In the model (1.3), we take the 


weight of the penalization term to be compactly supported in a ball, smooth, radial, repulsive 
and normalized to have mass 1. Specifically, we choose the function V : 1 -> 1 so that for 

V(s) :=nM) 

(Al) suppV = 6(0, r) for some 0 < r < 1, 

(A2) V € C' 1 (T n ), 

(A3) W(|x|) ^ 0 for all x £ T n , and 
(A4) f Tn V(x) dx = 1. 

Also, we take 

V : = 


a 

N 


for some a > 0. With these choices the energy is 0(1). Moreover, we assume that u p = 1, 
i.e., that the nanoparticles completely prefer the phase u = 1. With the above choices, the 


free energy (1.3) of a nanoparticle-polymer blend with IV-many round particles centered at 
x £ T" x ■ ■ ■ x F with radius r > 0 is given by 


(1.4) 


e,a,r,N(u;X.) = [ \Vu\ 2 dx + [ (u 2 - if dx 

o Jjn loe Jjn 


N r 

+ * § Lm 


V(\x — Xi\)(u — if dx. 


We consider the energy E e .a,r,N over functions u £ H 1 ( T n ) with f Jn u(x) dx = m and a set of 
fixed points x £ (T”)^ as centers of nanoparticles. Note that the constants in front of the first 
two terms in (1.4) are chosen so that these two terms together T-converge to the perimeter of 
the phase u = 1. namely to the first term of as e —X 0 in the L 1 (T")-topology (cf. |50]). 

1.4. Outline. Our main result in Section [2] states that for any absolutely continuous prob¬ 
ability measure ^ £ £P ac (T n ) we can find W-many points so that an appropriate extension 
of the energy E e)(Tir£> Ar e to L 1 (T n ) T-converges to the energy E^^ given by (1.1) as the num¬ 
ber of particles N e tends to infinity and the radius r e of each particle goes to zero as e —> 0 


(Proposition 2.1). To prove this result we exploit the above mentioned T-convergence of the 


periodic Cahn-Hilliard energy to the isoperimetric energy as the interfacial thickness e goes to 
zero (cf. 113 [,[50]) while approximating the measure /i by measures where the density is given 
as the weight in the penalization term of A classical consequence of T-convergence is 
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that a sequence of minimizers of the energies E e ,a,r e ,N e converges to a minimizer of the energy 
Ein ^ 1 (T n ) as e —> 0 (Proposition 2.3). Note that the energies E^g and E e;Cr admit 
global minimizers by the direct method of the calculus of variations for any /i G ^ ac (T n ), 
a > 0 and e > 0. This T-convergence result also extends to the copolymer model with the 
inclusion of the nonlocal repulsive interaction term (Remark |2.2[ ). 

In Section [3j we state regularity and criticality properties of local minimizers of the energy 
E up. Indeed, we prove that the phase boundaries of L l -local minimizers of E MjCr are regular 
provided the density of the measure fi is in L°° (Proposition 3.1). Moreover, under additional 
smoothness assumptions on the density (namely when the density is C 1 ) we present the first 
variation of the energy E^ i<r giving a necessary condition of criticality (Proposition |3.3[ ). 

The periodic isoperimetric problem has been the focus of much attention in the past (see 
e.g. [27-29,41 and references therein), and it is well-known that solutions of the isoperimetric 
problem possess phase boundaries of constant mean curvature. This, of course, may not 
be the case for the minimizers of the energy E^g. However, exploiting the regularity and 
criticality results of Section[3j in Section[4]we provide an example in two dimensions to illustrate 
how the minim izing pattern may be affected by the presence of the penalization measure fi 
(Example 4.1). Chosing dfi = ^xs(o,r) dx with appropriately chosen radius r (compared to 
the mass constraint m), we show that the lamellar pattern, which is the global minimizer of 
the classical isoperimetric problem [29| , ceases to be a minimizer of energy E^g for any a > 0 
with this choice of penalization. Moreover we discuss possible global minimizers of E^ iCr for this 
particular /r and prove that for all sufficiently large a, the global minimizer of E^g is given by a 
disk inside the support of fi (Propositions 4.2 and 4.3). These discussions also emphasize how 
the penalization affects the morphology (and geometry) of the phases (and their boundaries) 
(see Figure [T]). 

We conclude (in Section 5) by pointing out several directions for possible future studies. 





(c) a > 0 large 


Figure 1. Illustration of transition from a single striped lamellar pattern to 
a disk as a grows. Patterns (a) and (c) are the global minimizers of E^g for 
a = 0 and a > ao, respectively. Pattern (b) is a likely candidate for a minimizer 
with a > 0 but small. The gray disk represents the penalization region R(0, r). 


2. Large Number of Asymptotically Small Particles 

In this section we prove that given an absolutely continuous probability measure \i G 
<5P ac (T n ) the energy E Mi g appears as the asymptotic limit of the energy E^g^jy as the number 
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of particles approaches infinity, and the size of nanoparticles and the thickness of the phase 
transition go to zero simultaneously. Here the measure p gives the probability distribution 
of infinitely many asymptotically small nanoparticles in the e —> 0 limit. Indeed, we identify 


the energy (1.1) as the T-limit of the energy E e ,a,r,N in this asymptotic regime when r and 
N approach zero and infinity at certain rates, respectively, as functions of e. To this end, for 
given e > 0, let r e = r(e) and N e = N(e) be functions of e so that 


( 2 . 1 ) 


r e 


0 , 


N e 


Too and ^ r" < 1 


as e —> 0. 

For fixed Ay-many points xi ,..., xjv e GT" consider the energy E e a which extends the energy 
Ee,o-,r e ,Ay defined in (1.4) to L l (T n ) as follows 

f J T n \Vu\ 2 dx + ^ f Jn (u 2 - l ) 2 dx 
+ i 


( 2 . 2 ) 


Ee J o-('w) •— < 


N e r™ 


f B (x i ,r e ) V (\ x ~ x i\/ r z)( u - l f 


dx 


if u G H 1 (T n ) 

and J Tn udx = m, 


.Too 


otherwise. 


Note that here we rescale the weight as V £ = r~ n V(\x\/r e ) so that the assumption (A4) is 
satisfied for all r e —y 0 as e —> 0. 

Similarly we will extend the energy Ehowever, with an abuse of notation, we will still 
denote its extension to L 1 (T n ) by E^. Namely, for p G ^ > a , c (T n ) let 


(2.3) 


E/i,(T (vi) 


” \ fjn |Vu| T CT f Tn (u(x) - l ) 2 dp(x) 

< 


if u G BV(T n ),u = T1 a.e. 
and fj n udx = m, 


.Too 


otherwise. 


With these definitions we obtain that the family of energies E £]Cr T-converge to E^ in the 
topology of L 1 (T n ) as e —0. 


Proposition 2.1 (T-convergence of E £i(T ). Let p G ^ ac (T n ) be a probability measure that 
is absolutely continuous with respect to Lebesgue measure. Let the energies E £i(T and E M)(T be 
defined by \2.2\ and (2.3). Then there exists N e -many points xi,... ,xjy e in T n such that 

(i) (Lower bound) for any u G L 1 (T n ) and for any sequence {u £ } £ >o in L 1 ( T”) such that 
u e -r u in L 1 { T”) as e —> 0, we have 


lim inf E £ 

e ^0 


r(ti) ^ E Mi 0 -(u), 


and 


(ii) (Upper bound) for any u G LflT 71 ) there exists a sequence {w £ } £ >o in L l { T n ) satisfying 

v e —> u in L 1 (T n ) 


and 


lim sup E £!(T (u £ ) ^ E ^ a {u). 

e-S>0 


Proof. We will prove this proposition in three steps. 

Step 1. (Approximation). First, we will show that a measure p G ^ ac can be approximated 
in the weak-* topology by a collection of measures distributed on balls of radius r e > 0 with 
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weight V. Similar arguments appear in 49, Proposition 2.2] and [30, Lemma 7.5] for specific 
rates of convergence of r e and N e depending on the physical model. To begin, for p £ ^ ac (T n ) 
let p £ L 1 (T”) denote its density and suppose that 1/C < p(x ) < C for some C > 0 and for 
a.e. x £ T n . If not, we can consider measures p k defined via densities p k = (p A k) V 1/k and 
approximate p k by p%. Taking a diagonal sequence then will yield the result. 

Given e > 0 , define k e := 5 where [-J denotes the greatest integer less than its argument, 

and let I kf denote the index set { 1 ,..., 2 nfce }. Define the family of nested cubes {Qi}iei kc where 

(i) Qi C T" is a cube of side length 2 -Tlfc % 

(ii) Qi nQj = $ for i / j, 

(iii) T" = Uie 4 ^ or e > and 

(iv) for ei < e 2 , 4 C I k and {Qi} ie 4 C {Qijieiu- 


Also, define d e := 2 nfee , and along with ( 2 . 1 ), suppose that 

r e < N7 1/n < d e . 


For each i £ 4 f > let Nf = |_ N e p{Qi)\ . In each Qi select Nf points { x i-j}j=] that are almost 
equally distributed such that the distance between them is of order d e (Nf)~ 1 N. Noting that 
C/d £ ^ p{Qi) ^ Cd £ , this implies that for j ^ k we have 

(2.4) \x\j — x\ k \ > C » Cr e . 

Also, note that for each e > 0 the total number of points N e is 


Ne = E N i 


Let /x\ be given such that 

1 4 

d Vi( x ) : = ~k^2 V (\ X ~ x \j\/ r e)XB{xl j ,r t ){ x ) dx 


Nf 


and define 
(2.5) 


3 = 1 


:= w. S n 

iei he 


Clearly, p £ is a probability measure, as p\{ T n ) = Nf implies that p e ( TP 1 ) = 1. Moreover, by 
( 2~4]), for any * £ I kc , B(xf j ,r e ) C Qi and B(x e ip r e ) D B(x e ik , r e ) = 0 for j, k = 1,..., Nf. 

To prove that p e —» p as e —> 0 in the weak-* topology it suffices to show that 

lirnsup /j e (Q) ^ /i(fl) 
e^O 


for any closed set H C F. The weak-* convergence then follows from the Portmanteau 
Theorem (cf. [57j Theorem 1.3.4]). 
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For any closed set O C T", by the assumption (A4), we have that 
Me(^) = jr X 




1 


Nf 


€ 


N f rl 


XX/ V(\x-x\j\ /r e )dx 

iei ke j =l 

Nf 

X X / VQx-x^/r- e )dx 

iit. U J ^,re) 


= — y jv? 

iV £ ^ 

*£_4 e 

^£hc e 

onQf^0 


Let 5 > 0 be arbitrary. Define 

n e := J Qi. 

^-^fcg 

OnQ ^0 

Since the families {(5,;}j e / fc are nested, we have that D £l C fl £2 for e\ < 62 - Also, since D is 
closed, fl = He>o ^e; and f° r an y given 5 > 0 there exists €0 > 0 such that for all e < e 0 , 


X viQi) ^ M^) + <5- 


Therefore 

lirnsup f/, € (fi) ^ /r(D) + 5, 
e —>0 

and letting 6 —> 0 yields the result. 


Step 2. (Lower bound). After relabelling the points {x\-} found in Step 1, for any e > 0 we 
obtain a set of N e points Xi ,..., £ T n . For any u £ L 1 ( T n ) define the energy E eA given by 


(2.2) using the points {x',}X. 


and 


P e (u):=^-[ \Vu\ 2 dx+^-[ (u 2 -l) 2 dx, 

O Jjn loe JT n 


1 ^ r 

K e{u ) := —^ X / V’G® - ®<IAe)(“ - l ) 2 


so that E e,o(u) = P £ (it) + aK e (u). 






SHARP INTERFACE OF NANOPARTICLE-POLYMER BLEND MODEL 


9 


Let u € L 1 (T n ) and let {u e } e >o C L l (T n ) be a sequence such that u e —> u in L 1 (T") as 
e —> 0. We can assume, without loss of generality, that u = ± 1 a.e. Otherwise 


lim inf E c a {uA 

e —>0 ’ V ' 



l ) 2 dx = +oo, 


and the result of Part (i) follows trivially. Similarly, if Jj n udx ^ m , then for small e > 0, 
fjn u e dx / m and liminf^Q E t) a{u t ) = +oo. Therefore it suffices to consider only functions 
u 6 L 1 (T n ) satisfying u = ±1 a.e. and f Jn udx = m. 

Moreover, assume that —1 ^ u e sj 1 a.e. If not, we can consider the truncated functions 


u e 


-1 

on 

{x 

u e {x) < - 1 } 

U e 

on 

{x: 

- 1 < U t (x) 

1 

on 

{x: 

Ue{x) > 1 }. 


Clearly u e -A u in L^TP 1 ) as e —> 0. Also, P e {u e ) ^ Pe{u e ) and K e {u e ) ^ K e (u e ). Hence, we 
can use the truncated functions u e instead of u e . 


In [50| Section B], the author shows that liminfg^o P e (u e ) ^ P(u). Now we will show that 
a similar lower semi-continuity property also holds for the penalization term K e ■ Note that, 
by (2.5) we can write 


K e (u) = [ (u e - l ) 2 d/j, e (x). 

JT n 

Also note that since u e —> u in L 1 (T n ) and — 1 ^ u e ^ 1 we have that u e — > u in any L p (T n ) 
with p ^ 1. In particular, (u e — l) 2 —> (u — l) 2 a.e. as e — > 0. Then by the Egoroff’s Theorem, 
for any arbitrary 5 > 0 there exists a compact set G$ C T n such that (u € — l ) 2 — > (u — l ) 2 
uniformly on G$ and |T n \ < 5. Moreover, there exists eo > 0 such that for any e < eo 


(Ue - l ) 2 


(u 



6 

i^i 


for all ieGj, 

On the other hand, since u = ±1 a.e., (u — l ) 2 = 4 xa c where A = {x £ T n : u(x) = 1} and 
A c denotes the complement of A, namely, T n \ A. Since the set A has finite perimeter it holds 
that \dA\ = 0, and since p is absolutely continuous with respect to the Lebesgue measure we 
get that p(dA) = 0. Thus, A (and A c ), is a continuity set of the measure p. Since p e —>• P 
in the weak-* topology, this implies via the Portmanteau Theorem (again, cf. 

1.3.4]) that 


57, Theorem 


lim^ e (A c ) = p(A c ). 
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Combining these, we get that 

liminf / (u £ — l) 2 dp £ (x) ^ liminf / (u £ — l) 2 dji e (x) 

J T n e_5>0 Jgs 

= liminf f (^(u £ — l) 2 — {u — l) 2 + (u — l) 2 j dg, £ (x) 


^ liminf / (tt — l) 2 dg £ (x) — <5 


e ->0 


'G S 


= / (u — l) 2 dfi(x) — 5 
JG S 




(u — l) 2 dg,{x) — C5 


for some constant C > 0 independent of e > 0. Therefore letting 5 —> 0 yields 

liminf K £ (u £ ) ^ K(u ); 

hence, Part (i) follows. 

Step 3. (Upper bound). Let u £ L 1 (T n ) and let the points {xj}^! 1 C T n be given as above. 
Assume that u £ BV(T n ), u = ±1 a.e. and J T „ r(x) dx = m. Otherwise, = oo and by 

choosing v £ = u for all e > 0 the result of Part (ii) follows. Let the set A C T n be such that 


u(x) = 


1 if x £ A 

-1 if x £ A c . 


Let r = dA n dA c and assume that T £ C 2 . If not, one can approximate A by a sequence of 
open sets as in 50, Lemma 1] satisfying the condition that dA & is of class C 2 . Then 


one can use the sets dAk 0 dA c k instead of T to prove Part (ii) and then pass to a limit using 
a diagonal sequence. 

Define the signed distance function dr : T n —> M by 


dv{x) ■■ = 


dist(x, T) 

— dist(x, T) 

and the sequence of functions g £ : M —> M by 

-1 

- 2,/S) - 1 

9e(t) ■■= < 

t + 2-y/e) + 1 



if x £ A c 
if x £ A, 

if t > 2^/e 
if y/l ^ t ^ 2y/e 
if \t\ ^ sje 

if — 2 s/e ^ t ^ — yje 
if s < —2y/e 


where the function z(t) solves the ordinary differential equation 

^ z 2 — 1) subject to z(0) = 0. 

In [50] the author shows that the function defined by 
(2.6) v £ (x) := g e (dr(x)) + rj e , 
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where rj e is an additive constant that is of order O(e), is in H 1 ( T n ), and satisfies the mass 
constraint f Jn v £ (x) dx = m for all e > 0. Moreover, 

v e —>• u in L 1 (T n ) as e —> 0 and lirn sup (w e ) ^ P(u). 

e ->0 

We will use the family of functions u e ’s to prove that a similar limsup inequality also holds 
true for the penalization term K e . To this end let 

A e + := {x € T n : g e (dr(x)) = 1} and At := {x € T n : g e (dr(x)) = —1}, 

and let T e denote the transition layer 

T e :={xeT n : - l<g e {d r (x)) < 1 } 


Define 


Then, for any e > 0, T n = 
Using the fact that r] t = 
we get that 


r e — r 

L + ■— 1 < 


n A and T1 := T e n At 


A £ + U U TL U At, A = r*. U A e + and A c = Tt U At. 

= 0(e) and g € (dr(-)) G L°°(T n ), for the functions v e given by (2.6) 


Ke{Ve) = (Ve ~ l) 2 dfi e (x) 

J T n 

= J ^ (g e {dr{x)) - l) dg e (x) + J (g e (d r (x)) - lj dg e (x) + 4 J dg e (x) + 0(e) 

= (ge{dr(x)) - l) dg, e (x) + ^(g e (d r {x)) - l) 2 - 4^ dfi e (x) + 4^ dfi e (x) + 0 (e) 

^ 4// e (T^_) + 4 f d/j, e (x) + 0(e). 

Ja c 

Let 5 > e > 0 be fixed but arbitrary. Then C Tj_, and for T <5 = U Tl we have 

K e (v e ) ^ 4/i e (Ts) + 4 f dfi e (x) + 0(e). 

Ja c 

Again, by the Portmanteau Theorem the weak-* convergence of to g is equivalent to the 
fact that limsup^o /^(T^) ^ as T^ is closed in TL Moreover, since \dA\ = 0 and A 

is a continuity set for the measure /i by its absolute continuity with respect to the Lebesgue 
measure, as in Step 2, we have that lim e _ > o Ah(A c ) = g(A c ). Therefore, 


limsup Kfiyt) ^ Ag(Ts) + 4 / dg,(x) ^ C5 + K(u) 
e —>0 J A c 


for some constant C > 0 independent of e > 0. Letting i A 0 and combining this with 
limsup^Q P e (v e ) ^ P{u) we obtain the result of Part (ii). □ 
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Remark 2.2. (Nonlocal perturbations and the diblock copolymer model) Define the nonlocal 

perturbations of the functionals E e a and E^, a respectively by 

(2.7) 

/t« |Vu| 2 dx + f T „(u 2 ~ 1 ) 2 dx if u G H 1 ( T n ) 

Eil £ i / B ( Xi , re ) V(|x - Xi\/r e )(u - l) 2 dx and f Jn udx = m, 


Ee,<7, 7 (^) •— < 


+7 fj n Jir n G ( x ) 2/) (' u ( x ) “ m ) (“(2/) ~ m ) dxd V 


+oo 


otherwise, 


and 

( 2 . 8 ) 




5 It" + a fjn( u ( x ) — l) 2 dy{x) if u G i?D(T n ), u = ±1 a.e. 

+7 fjnfjn G(x, y) (u(x) - m) ( u(y) - m) dxdy and f Tn udx = m, 


+oo 


otherwise 


over functions u G L 1 (T"). Then a standard conclusion of T-convergence is that the conver¬ 
gence is stable under continuous perturbations, i.e., E e)(Ti7 T-converges to E M)(Ti7 as e —> 0 in the 
L 1 (T”)-topology (cf. [l7j Proposition 6.21]). 


Another classical consequence of T-convergence is that the limit of a convergent sequence 
of energy minimizers minimizes the limiting energy. The proof of the following proposition is 
quite standard and can be adapted, for example, from the proof of [50| , Theorem 1], 

Proposition 2.3 (Limit of a sequence of minimizers). Let y G ^ ac (T n ) and let x \,..., xjv £ be 
such that the sequence of measures {y £ } e >o C ^ > ac (T n ) defined via the densities 

1 Ne 

( 2 - 9 ) Xi \/ r JXB( Xi ,r e )( x ) 

e 6 i= 1 


converges to y in the weak-* topology of £P(T n ) as e —> 0 and N e —> + 00 . Suppose u € —>• u in 
L 1 (T n ) as e — > 0 where, for any e > 0 , u e minimizes the energy E e>(T for any a > 0. Then u 
minimizes the energy E /t , CT over BV{ T") with u = ±1 a.e. and J Tn udx = m. 


Remark 2.4. Note that given any measure y G ^ ac we can find lV e -many points as in the Step 
1 of the proof of Proposition |2.1| so that the probability measures y e defined via the densities 


(2.9) converge to y in the weak-* topology. 


Remark 2.5 (Compactness of a sequence of minimizers). Using the polynomial growth of the 
double-well potential {u 2 — l) 2 and the compactness of L>D-functions in L 1 


(cf. 


24 


Theorem 

1.19]) we can easily conclude that if {u £ } e> 0 is a sequence of minimizers of the energies E fj0 - 
then there exists a subsequence {u f . such that u ek uq in L l (T n ) as ej —> 0 for some 
function uq G L 1 ( 


(cf. 50, Proposition 3]). 
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3. Properties of Local Minimizers of E 


/i,(T 


Independent of its connection to nanoparticle-polymer models, the energy E ^ also piques 
one’s interest as a rather simple extension of the classical periodic isoperimetric problem where 
one tries to minimize the perimeter of a set of fixed mass with respect to a penalization term 
determined by a fixed probability measure. Indeed, as a model for pattern formation, the 


minimization of (1.1) sets up a basic competition between short-range effects of the perimeter 
term and possibly long-range penalization via the choice of the measure p. The interplay 
between these competing terms appears in properties of local minimizers such as regularity, 
criticality and stability. 

For /i E ^ac(T n ) let us denote its density by p E L 1 (T n ) for the remainder of this section, 
i.e., let 

dfi(x) = p(x) dx. 

With a slight abuse of notation we will also denote by E // fJ the energy defined on sets of finite 
perimeter. Namely, we consider the problem 


(3.1) 


locally minimize E^ a {A) : = 


I VxaI +4 a p[x) dx 

J M 


over sets of finite perimeter A C T” such that |A| = (1 + m )/2 for any given m E (—1,1), 
a > 0 and p E L 1 (T n ). Note that this formulation is equivalent to minimizing (1.1) subject to 


the constraint f T „ u(x) dx = mn. Let us also note that a set of finite perimeter f i is an L 1 -local 


minimizer of E^ a if 


(3.2) 

for some 5 > 0. 


E^ l(T (I2) < E ^(A) provided 


Ixn - Xa\ dx < 6 


As we noted in the introduction the minimization problem (3.1) is in the spirit reminiscent 


of finding minimal boundaries with respect to an obstacle set. These obstacle problems tackle 
the following minimization problem: 


(3.3) 


minimize 


IVxaI 


over sets of finite perimeter AcF such that L C A where the obstacle L is a given fixed set 
of finite perimeter. Note that here the admissible sets are not mass constrained as they are 


in our problem. We believe that such an isoperimetric obstacle problem is equivalent to (3.1) 
in the limit a —> oo when |A| = \L\; however, when minimizing E^ )CT (A) we do not explicitly 
restrict the admissible patterns A to those which must contain supp p. 


We note that the regularity of phase boundaries of the obstacle problem (3.3) were estab¬ 
lished in 54, Section 3] depending on the boundary regularity of the obstacle set L. For our 
problem (3.1), on the other hand, we prove that the regularity of the phase boundaries are 


determined by controlling the L°°-bound of the density p rather than the smoothness of the 


boundary of its support. The issue here is to control the excess-like quantity (3.10) that mea¬ 


sures how far a set is from minimizing perimeter in a ball in terms of the radius of that ball. 
Indeed, we show that if p has bounded density p then for any a > 0 the penalization term 
can locally be controlled by the perimeter term and we can conclude by the well-established 
regularity theory for the isoperimetric problem that the phase boundary of a local minimizer 
of E^ )(r (A) is of class C 1,a . The essential elements of the regularity result are already contained 
in the works of others in the similar context; however, we are unaware of a particular result 
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that applies to our setting specifically. Hence, for completeness, we present here a proof for 
the regularity of phase boundaries. 

Proposition 3.1 (Regularity of Phase Boundaries). If p £ L°°(T n ) and Ll C T n is an L 1 -local 
minimizer of (3.1), then d*Ll is of class C 1,a for some a £ (0,1) and R s (dLl \ d*Ll) = 0 for 
every s > n — 8 where d*Ll denotes the reduced boundary of LI and'H s denotes the s-dimensional 
Hausdorff measure. 

The proof of the regularity of local minimizers of E^(H) rely on the following technical 
lemma proof of which can be found in 23, Lemma 2.1]. 

Lemma 3.2. Let L be a Borel set, and let D be an open domain such that f D |V X i,| > 0. 
Then there exists positive constants ko and Iq depending only on D and DHL such that for 
all k with \k\ < k$ there exists a set F such that F = L outside of D and 

|T| = \L\ + k 


Id 


|Vxf| ^ / \V X L\ + lo\k\ 


ID 


[ \xf ~ Xl\ dx ^ lo\k\ [ |V XL |. 
J D J D 


Now we prove the regularity result proceeding as in [51} Proposition 2.1] 


Proof of Proposition 3.1 Let LI be an L 1 -local minimizer of (3.1) and xq £ d ll be arbitrary. 
Let D CC T" be such that xq 0 D and f D |V X n| > 0. For L = LI in Lemma 3.2 there exist 
two constants k o and Iq that depend only on D and find. Using these constants fix R > 0 
such that 


(3.4) 


OJnR F ko, 


1 + fo / |V X n| 

Jd 


oj n R n < 5 and Br(xq) flfi = 


where co n is the measure of the unit n-ball and <5 is as in (3.2). 

Let minimize the perimeter in Br(x q) subject to the boundary values of fi, i.e., 


’ B r {x q) 


l V Xol ^ 


'Br(xo) 


|V X a| 


for all A Cl" such that A \ Br(xq) = Ll\ Br (x q) . 


Since fifl D = Ll H D the result of Lemma 3.2 holds true with the same constants k$ and 


Iq if we replace LI by LI. Hence, for k := |H| — |H| ^ ai n R n < ko by the choice of R > 0, there 
exists a set G such that G = LI outside D and 


(3.5) 

(3.6) 

(3.7) 


/ I XG 
J T™ 


id 


\G\ = |fi| = m, 

|V XG | ^ / |Vx„| +CR n , 

JD 


- xd dx ^C 0 R n <5, 


where (3.7) follows from (3.4) with Co := (l + lo J D |V X n|) oo n . 

By (3.5) and (3.7), the set G is an admissible competitor for the energy E ^^(A)-, hence, 

/ |V X n| + 4 ct / p(x)dx ^ / |V XG | + 4<r / p(x)dx. 

J T n Jn c J T n Jg c 
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Noting that Cl \ Br(xq) = Cl\ Br(xq) and G\D = Cl\D, and using (3.6) we get that 

/ 

JT”\(D\JB r (xo)) 


|Vxn| + / !Vx^l+ / |Vxn|+4a [ p(x)dx 

JD J Br(xq) J £l c 


ID 

< / |Vxg| + / |VxgI 

JT"\(DUBjj(io)) JD 


+ / IVxg| + 4cr / p(x)dx 

Jb r (x o) 4g c 

/ |Vxn| + / IVxgI 

/p\(DUB fl (io)) 4G 

+ / |VXnl+ 4cr / P(z)cte 

J Br(x o) jG c 

/ iVxnl + / IVxgI 

+ f |VxqI+4o- f p(x) dx + C R n 
J Br(xq) JG c 


€ 


for some constant C > 0. Hence, we have 
(3.8) 


(3.9) 


/ |Vxn| — / IVXqI ^ 4cr f/ p{x)dx- [ p(x) dx\ + CR n . 

'B r (x o) ^Sfl(xo) V4 g c 4^ / 

On the other hand, using (3.7) and the fact that p £ L°°(T n ) we obtain 
p(x) dx — p(x) dx = I p{x) dx — p(x) dx 

JQ. c Jn Jg 

= / P(x)(xn ~ Xa){x)dx 
J T n 

^ l|p||L°°(T")IIXf2 - XgIIl^t-) ^ 072". 


/G c 


Combining (3.8) and (3.9), we get that 

IVxnl-/ |Vxnl<Oi 2 ". 

J Bntxn) 


(3.10) 


'B r {x o) 


Property (3.10) states that the boundary of the set Cl is almost area-minimizing in any ball. 


With this property, the classical regularity results of 40,54] apply, and we can conclude that 
d*Cl is of class C 1,a , with B s {dCl \ d*Cl ) = 0 for every s > n — 8. □ 


With the regularity of phase boundaries at hand, under further smoothness assumptions on 
the density p we have the following necessary condition of local minimality. 

Proposition 3.3 (Criticality Condition). If p £ C' 1 (T") and u is an L 1 -local minimizer of the 
energy E Mi(T , then 

(3.11) (n — 1 )H(x) — Aap(x) = A for all x G dA 

for some constant A where H : dA —> M denotes the mean curvature of dA and A = {x: u(x) = 
1} as before. 
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Proof. Suppose p £ C 1 {T n ) and let u be an L 1 -local minimizer of ELet A = {x\ u(x) = 1}, 
and let ( £ C°°(dA) such that f QA ((x) d'H n ~ 1 (x ) = 0. 

To compute the first variation of the energy E^ iCT , we view it as a set functional given by 


Let 

X(x) = ((x)v(x) on dA 


(3.1) and proceed as in 14,52 . 


where v denotes the outer unit normal to dA. Then, clearly, 


(3.12) 



X-vdU n ~ x {x ) = 0. 


Let T : T n x (—r, r) —> T n solve 
(3.13) 


T(x, 0) = x, 


for some r > 0. Define 

A t := 

Invoking Proposition [3Tj we easily see that for the family of sets {4 t } tg (_ r r ) for some r > 0 
we have that 


XA t 


dA t is of class C 1 , and 
-A xa as t -A 0 in L l (Y n ). 


Moreover, 


(3.14) = I + tVX + ]a 2 XZ + o{t 2 ), 

where Z := d 2 ^/dt 2 \ t= o is given with z-th component = Y^j=\ Xxjx^K and 


(3.15) 


d_ 

Ft 


JT = traceVX = divX, 
t =o 


where JT denotes the Jacobian of T. 
get that 


Hence, using (3.12) and the Divergence Theorem, we 


d_ 

dt 



d_ 

dt 


t= o. 


JT dx = 0, 


i.e., the family of sets {At}t£(- T .r) preserves the volume of A to first order. Therefore this 
family of sets is an admissible class of perturbations of A to compute the first variation of 

E ^(A). 

Define the functions U(x,t ) by 


(3.16) 


U(x, t) 


1 if x £ At, 
— l if x £ A^, 


and note that a function u £ BV{T n , {±1}) is said to be a critical point of the energy E fli(T if 
d/dt| t= oE MiCr ([/(•, t)) = 0 for every U(x,t ) defined via the admissible family {A t }. 
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Consider the energy 

E nAU(;t)) = \\ |W(-, t)\ + a [ (U(x, t ) - l) 2 p(x) dx 

Z Jjn Jjn 

= : P{t)+aK{t). 

In 114], the authors show that 

(3.17) P'(0) = (n — 1) [ H(x)C(x)d'H n ~ 1 {x) 

JdA 

where H denotes the mean curvature of dA. Now we are going to compute iC(0). Note that, 


by (3.13), 


K(t) = 4 / p{x) dx 
J A? 


= 4 


= 4 


/ p(x) dx — 4 p{x) dx 
J T n J A t 

/ p(x) dx — 4 / p(4/(x,t))J4’(x,t)dx 
J T n J A 


Therefore, since p G C 1 (T”), 

/• O r\ 

(3.18) K'(t ) = -4 J Vp(^(x,t)) — 4>(x,t)J4>(x,t) + f(V(x,t)) — {J4>(x,t)) dx. 


Hence, by (3.13), (3.14) and (3.15), using the Divergence Theorem we get that 

K'( 0) = —4 [ X7p(x) ■ X(x) + p(x) div X(x) dx 


=—4 di v(p(x)X(x)) dx 

Ja 


=—4 p(x)(X(x) ■ v{x)) dK n (x). 

JdA 


Combining this with (3.17) we get that 

[(n — 1 )H{x) — 4ap{x]\C,{x) d'H n ^ 1 (x) = 0, 


(3.19) 


IdA 


i.e., there exists a constant A such that 

(3.20) (n — 1 )H(x) — 4 crp(x) = A 

for all x G dA. 


□ 


Remark 3.4. Proposition 3.3 holds locally in case p is piecewise C 1 . That is, if we assume p 
is C 1 except on a smooth submanifold (on which it or its derivative is allowed to jump), the 
curvature condition (3.11) holds at regular points of p. This observation follows by noting that 
the weak form (3.19) continues to hold for ( supported in each component of the set of regular 
points of p, and that by appropriate choices of ( we may conclude that the Lagrange multiplier 
A is independent of the component. 
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Remark 3.5. Since the boundary of A = {x: u(x) = 1} of an L 1 -local minimizer u of the energy 




is of class C 1,a by Proposition 3.1, we can express the reduced boundary d* A locally as 


the graph of a C 1,a function ip on a ball B C T n 1 . Then, the first variation (3.20) implies 
that 

V<p{x') 


(n — 1) V • 


= 4ap(x', f(x')) + A for x' E B. 


V 1+ l v ^( x ')l 2 / 

As the right-hand side is of class C 1 , by standard elliptic regularity we obtain that ip E C 3,a . 
Hence, the boundary of A = {x: u(x) = 1} is of class C 3,a for some a > 0. 

Remark 3.6. Note that the condition ( 3.1 1[ ) is a sufficient condition for u to be a critical point 
of the energy E Mj(T with respect to L 1 -perturbations. 

Remark 3.7 (Second Variation). If we further assume that p E C 2 ( T"), then a necessary 
condition for L 1 -local minimimality of u is given via the second variation of the energy E ^ 
around the critical point. Namely, 


(3.21) 


(|v 9 aCI 2 - II^9a|| 2 C 2 


dW 1 - 1 - 4a 


(yp-v)Cdu n -\x) ^ o 


for any smooth ( : dA —> M satisfying J dA ( dH n ~ 1 = 0. Here VsaC denotes the gradient of C 
relative to the manifold dA , Bqa denotes the second fundamental form of dA and v denotes 
the unit normal to dA pointing out of A. The computation of (3.21) follows by adapting the 
calculations in 14, Theorem 2.6]. 


In the absence of nanoparticles (when a = 0) an important result regarding the local rnini- 
rnizers of the nonlocal isoperimetric problem related to the energy (2.8) is given in [I . Here the 
authors prove that strict stability in the sense of positive definite second variation of critical 
sets is a sufficient condition of isolated local minimality with respect to the L 1 -topology. We 
believe that the techniques introduced in 11] can be adapted for the functional E M](T to conclude 
that strict positivity of (3.21) implies local minimality in L 1 . 


4. An Example in Two Dimensions 

Depending on the distribution of nanoparticles and the strength of penalization via a one 
can modify the phase morphology of block copolymers and effectively prescribe the location 
and shape of the phase transitions. Even with a given measure p E £A-&c^ n ') describing the 
particle distribution there are many possible critical patterns for the energy E MjCr , depending 
on the strength coefficient a > 0. Indeed, the penalization term can act as an attractive or 
repulsive term depending on the choice of p via its pinning-like quality. The rigidity of the 


results in Propositions 3.1 and 3.3, on the other hand, limits the possibilities for critical and 
minimizing patterns. In this section we will provide such an example in two dimensions, i.e., 
on the 2-flat torus T 2 . Exploiting these rigidities, the example below shows that for a certain 
choice of p and when the mass constraint m is restricted to a certain range, for any <r > 0 
the global minimizer of the energy E^^ is geometrically quite different than the solution of the 
isoperimetric problem, i.e., when <7 = 0. 

Example 4.1. For T 2 = [—1/2,1,2) x [—1/2,1/2) with periodic boundary conditions, any 
m E [0,1 — 2/7r) , and any fixed r > yj (27r)” 1 (l + m) let p E A^ ac (T 2 ) be defined via the 
density function 


PW ■= ^2 XB(0,r)(x) 
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for X G T 2 . 

By the direct method in the calculus of variations, there exists a global minimizer of for 
any a > 0 and for the measure p, defined as above. Let also A denote the set {x G T 2 : uq{x) = 
1 }. Since f j2 uo(x) dx = m, we have that |A| = (1 + m)/2 and \ A C \ = (1 — m)/2. The problem 
thus reduces to find a set AcT“ with area |A| = (1 + m)/2 which minimizes 

E M , ff (A) = Per x 2 (A) + 4a (l - ^ r)l ) . 

That is, A should have as small a perimeter as possible, while maximizing its intersection with 
the nanoparticle domain B(0,r). We note that for m G [0,1 — 2/ 7 r), (1 — m)/2 < (1 + m)/2 
and by the choice of r as above, we have that |B(0,r)| > (1 + m)/ 2; hence, An B(0,r) 0 for 

any admissible set A. 

When a = 0 the material is nanoparticle-free, and the energy reduces to the classical isoperi- 
metric problem on T 2 . For our choice of mass m G [0, 1 — 2/ 7 r), the unique minimizing config¬ 
uration with <7 = 0 (up to translation) is the single striped lamellar pattern 


ul(xi,x 2 ) = 


if x\ G ( 
if X! <£ ( 


— 1 —m 

4 ’ 

— 1—m 


l+m \ 

4 ) ’ 
1+rnA 
4 ) > 


— 1—m L+m'i 


)}. (See Figure 1(a) in the Introduction.) As 


with associated set Aj j : = {aq G ( 
noted above, by the choice of radius r, any translate of Al must intersect the nanoparticle site 
.B( 0 ,r); as we will see below, this will imply that the lamellar pattern cannot be the energy 
minimizer for any a > 0, and in fact it will no longer be critical for the energy E Mj(T . 


The shape of the minimizer is constrained by the curvature equations (3.11) which are 
satisfied by any critical c onfi guration. Indeed, since the penalization density p is uniformly 

we conclude that dA is of class C 1,a for some a > 0. On the 


3.1 


bounded, by Proposition 
other hand, as p is constant on each of int(An B(0, r)) and int(Ani? c (0, r)), where int denotes 


the interior of these sets, it is trivially differentiable, and the formula (3.11) is locally valid (see 
Remark |3.4| ) Thus 

H(x) = A for x € dA n int(L> c (0, r)), and 


(4.1) 


H(x) = + A 


7rr^ 


for x G dA fl int(R(0, r)) 


for some constant A. We note that H(x) denotes the signed curvature, and it is piecewise 
constant. In particular, when H > 0 the domain A lies inside a circle of radius 1 /H, while for 
H < 0, A is exterior to a circle of radius 1/\H\. 

We may immediately confirm the claim made above, that the lamellar configurations, con¬ 
sisting of translations of ul , cannot be minimizers (or even critical) for any a > 0. Indeed, by 
our choices of parameters m,r, any translation of Al intersects intL>( 0 ,r), so the curvature 
condition is violated in the intersection. 

Another observation which follows directly from the curvature conditions (|4.1|) is that any 
ball An = B(p,R), with R := yj(2 n) 


stationary for (See Figure 1(c) 


1 (l + m) and p G T 2 such that A ]j C B(0,r), is 
below.) This configuration is a local minimizer for any 


( 7 ^ 0 , and in fact it is the global minimizer for all sufficiently large a: 


Proposition 4.2. There exists uq = ao(m,r) such that for all a > or, Ar (defined above) is 
a global minimizer of E MjCr . 
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We defer the proof of Proposition 4.2 to the end of the section. 

The question is then what is the geometry of minimizers for small positive values of a. As the 
energy depends continuously on the parameter a > 0 , for small values of a > 0 we expect that 
the global minimizer of the energy E MiCT is L 1 -close to a lamellar pattern when 0 ^m<l — 2 / 7 T 
(Figure l(b)| ). Below we propose a possible geometry for minimizers for small a. Although we 
cannot describe them completely, a minimizer which is not a disk must be “stripe-like” in the 
sense that it must exploit the topology of T 2 : 


Proposition 4.3. Let m £ [0,1 — 2 / 7 r] and r > y 2 (2tt) 1 (1 + m), and A C T 2 corresponding 
to a minimizer of E //iCT . Then: 

(i) If A is contractible in T 2 , then A is a ball of radius R = \J ( 27 t ) _1 (1 + m) with 
A C B(0,r). 

(ii) A c cannot be contractible in T 2 . 

Proof. First assume A C T 2 is contractible. We lift T 2 to M 2 , its universal cover. Contractibil- 
ity in the torus implies that the lifting of A consists of a periodic array of disjoint compact 
components A C M 2 , each with area |A| = (1 + m)/2. By the classical isoperimetric inequal¬ 
ity, each component has perimeter Per T 2(<9A) ^ Per T 2 (i?#), with equality if and only if the 
components are disks of radius R. By placing a periodic array of disks of radius R inside the 
array of translates of the nanoparticle site B( 0, r), we obtain a configuration which has smaller 
perimeter and which optimizes the penalization term, and thus has smaller energy than A, 
unless A were also a disk of radius R contained in B(0,r). Thus (i) is verified. 

The case of AT C T 2 contractible is similar. Since A c lifts to a periodic array of compact 
sets in M 2 , and Per T 2(cL4 c ) = Per T 2 (dA), we may conclude that a disk of radius R again has 
smaller perimeter than A. By locating the disk inside B(0,r) the penalization term in E AIj(T is 
optimized, so again the disk has strictly smaller energy than any domain with A c contractible. 
This proves (ii). □ 


Using Proposition 4.3 and the curvature condition (4.1) we may illustrate some configura¬ 


tions which are candidates for the minimizer, and eliminate certain others. Supposing that the 
minimizer is not a disk inside B(0, r), we may assume that both A and AT are not contractible, 
and hence each intersects both B(0,r) and B c (0,r). By the criticality conditions (4.1) we see 
that dA has to be a union of arcs of circles and straight lines as its connected components 
have constant curvature in two dimensions. Also, note that the curvature of dA inside the ball 
.B(0,r) has to be greater than the curvature of dA on £> c (0,r). Thus, dA does not consist of 
a union of straight lines inside B(0,r) and arcs of circles outside of B(0,r). 

Since dA is of class C 1,a constant curvature components of dA meet tangentially on dB{ 0, r). 
Therefore, dA can not consist of a union of an arc of a circle inside B(0, r) connecting to another 
arc of a positively curved circle outside of B{ 0, r), since two points and tangents at those points 
determine a circle uniquely, and two circles with different positive curvatures cannot meet at 


two points tangentially. Therefore the Lagrange multiplier A in (4.1) cannot be strictly positive 


since for A > 0 the components of dA would consist of positively curved arcs of circles which 
is not possible. Therefore we may assume that either A = 0 or A < 0. 


Band aid patterns. Suppose first that the Lagrange multiplier A = 0 in (4.1). In this case 
the domain consists of arcs of circles inside B(0,r) and straight lines outside of B(0,r). Note 
that, by periodicity of the domain, the straight components of dA in B c ( 0, r) must be parallel, 
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and hence they must meet 95(0, r) at semicircles inside 5(0, r). Such patterns we will refer 
to as band aid patterns (see Figure [ 2 ]). 



(a) One band aid 



(b) Two band aids 



(c) Slant band aid 


Figure 2. Band aid patterns are stationary, but are not global minimizers for 
m £ [0,1 — 2/7r) for any a > 0. Here again the gray disk depicts the penalization 
region 5(0, r). 


By adjusting the radii of the semicircles, we may match the area constraint and so these 
do represent stationary points of t he e nergy E Mi(T . However, as the bandaid patterns are all 
contractible in T 2 , by Proposition 4.3 they can not be minimizers for any a > 0, and thus 
A = 0 is not achievable for a minimizer. 


Concave/convex strips. Suppose the Lagrange multiplier A < 0 in (4.1). Then A lies inside 
of arcs of circles of radius R 2 inside of 5(0, r), and outside of circular arcs with radius Ri 
outside of 5(0, r). (See figure [3j) Moreover, the curvature condition (4.1) relates the radii to 
the parameter a via 


(4.2) 


2(7 _ 1 1 

7 rr 2 R\ R 2 


Lemma 4.4. If A is a minimizer, then R\ > ' — r. 


Proof. To verify the lemma, assume instead that Ri ^ \ ~ r - The connected components of 
A c \5(0, r) are then either contained inside circles of radius R\ which are disjoint from 5(0, r), 
or are bounded by disjoint arcs of this radius which connect to 5(0, r) at two points on 95(0, r). 
In the latter case, we note that by lifting to M 2 , the distance between adjacent images of the 
nanoparticle domain 5(0, r) is 2(^ — r) > 2R\. Thus, the images in M 2 of the components 
of A c are compact, and hence A c is contractible in T 2 . This contradicts Proposition 4.3, so 
therefore the lemma must hold true. □ 


We may now present our guess for stripe-like minimizers when a > 0 but small. We choose 
the inside radius R 2 so as to create a concave/convex stripe pattern as in Figure |3j In order 
to do this, it is necessary that the inside radius R 2 > r; otherwise, the circular arcs inside 
5(0, r) may not connect across the nanoparticle zone, and a curvy band aid pattern would 
result (see figure [ 4 }) As such a pattern is contractible in T 2 , it cannot be a minimizer. The 
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Figure 3. The set A is enclosed by arcs of circles where dA is negatively 
curved outside of B{ 0, r ) and positively curved inside -£>(0, r). The penalization 
region B{ 0, r) is highlighted by the gray disk. 


concave/convex stripe thus requires a lower bound on both 
realized for 


7rr 2 f 1 1 


7r r 


and hence can only be 


a = 


+ 


< 


2(1 - 2r) 


i?l i?2_ 

The exact values of R\ , i ?2 (and the centers of the constructing circles) will be also determined 
by the area constraint |H| = (m + l)/2 and the requirement that the resulting curve is C l ' a . 


Remark 4.5. We conjecture that when a > 2 ( 1 —2r) th en the minimizer must be a disk of radius 
R, inside the nanoparticle region B( 0, r). The variety of concave/convex regions which may be 
drawn is great (and is not restricted to shapes depicted in Figures [3] and |4j) so the optimum 
value of (To in Proposition |4.2| remains an open question. However, we observe that as <7 gets 


larger, the radii R 2 of arcs within B(0, r) must get smaller in order to satisfy (4.2) (given 
Lemma 4.4), and hence the area contained in the nanoparticle region is eventually insufficient 
to reduce the penalization term in the energy. (See Figure[4j) This observation forms the basis 


for our proof of Proposition 4.2 


We conclude with the proof that the disk of radius R gives the global minimizer for suffi¬ 
ciently large < 7 . 


Proof of Proposition^//^ Let A be the set associated to a global minimizer of E^, and set 
j3 := R 2 /Ar 2 < 1/4. If A is not a disk of radius R, then by Proposition 43 ifl B(0,r) c 7 ^ 0, 
and so dA consists of arcs of circles of radius R\ > 0 outside B(0,r), and of radius R 2 > 0 
inside B(0,r), satisfying (4.2). By Lemma 4.4, we have 


1 2cr 2 

i ?2 Tr 2 1 — 2 r ’ 


and thus there exists a\ = <7i(r) so that for all <7 > <j\ we have R 2 < fir/ 2. 

We now claim that for all <7 > < 71 , Ad B(0,r) lies inside a disjoint collection of circular arcs, 
each of which lies within distance fir of dB(0,r). Indeed, dAd B(0,r) consists of circular arcs 
of radius R 2 < fir/2 (by the above estimate,) either connected to dB(0,r) at the endpoints, 
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(a) Convex/concave pattern 
with contractible 
boundary 



(b) Convex/concave pattern 
with uncontractible 
boundary 


Figure 4. The sets A shown are enclosed by arcs of circles where dA is nega¬ 
tively curved outside of B(0, r) and positively curved inside B(0,r), with small 
radius R 2 < r. The left hand set yields a contractible region which cannot 
be a minimizer. The image on the right is unlikely to be a minimizer, as the 
penalization term will be large for a > 0. 

or as disks of radius R 2 contained in the interior of B(0,r). The arcs which contact dB(0,r) 
lie within distance f3r of dB(0,r ) by the bound on R 2 , so it remains to consider interior disks. 
First, assume that several such disks are contained in the interior of B(0,r). By the classical 
isoperimetric inequality, the perimeter of A would be reduced by replacing these by a single 
disk with the same total area, with no change to the penalization term, and thus reducing the 
total energy. However, this contradicts the minimality of A, and thus there can only be a single 
disk of radius R 2 inside B( 0, r). By translating this single disk to be tangent to dB( 0, r), the 
energy of A remains the same, so we obtain a minimizer with all components of A n B(0,r) 
within distance (3r of dB(0,r), as claimed. 

By the claim, A FI B( 0, r ) lies within an annular region B( 0, r) \ B( 0, (1 — j3)r) of thickness 
f3r. In particular, 

| A n B( 0, r)| < 7T (r 2 — [(1 — f3)r] 2 ) < 2irr 2 /3. 

We may then compare the energy of A to that of the single disk Ar C B(0,r) of radius R = 
y/ (27t) _1 (1 + m). By the isoperimetric inequality on T 2 , Per T 2 (H) ^ Per T 2 (Hi) = 2 < -kR 2 , 
given our choice of parameters. Thus, 

E Mi(7 (H) - E ^(Ar) >2-2KR + Aa(%- 

\r z 7 r 

/ d2 

^ - [2-kR -2]+4a (-J -2/3 

c>2 

= - [2irR - 2] + 2a— 

> 0 , 

2 

for all er > max{cxi, ttR — 1)} := ao(m,r). Thus, for all a > <ro, the minimizer must be a 
disk contained inside B(0, r). □ 
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Remark 4.6 (Small translations). Note that unlike the function ul, the function ur is unique 
only up to small translations, i.e., any translate ur )CL of ur defined by 

. . f 1 if x G B(a, R), 

1—1 if x f. B{a , R), 

for any a £ T 2 with |a| < r — R, we have that E^^ipiR^) = E A j i 0 -(rr) if R < r. The energy 
E^o-, though, is not translational invariant in general. This also reflects the “pinning” effect 
of the penalizing measure p. 


Remark 4.7. We believe that the results of Example |4.1| (Propositions |4.2| and 4.3) can be gen¬ 
eralized easily to the case when the penalization measure is given by an indicator function p sat¬ 
isfying (i) J T „ p(x ) dx = 1, (ii) B(p, R) CC supp p for some p € T n where R = \f (27r) -1 (l + m), 
and (iii) | suppp| > (1 + m)/ 2 . 


Remark 4.8 (The effect of a). The effect of the penalization term in E /(i(7 is more rigid than the 
effect of the nonlocal perturbation in E^^-y given by (2.8). Indeed, in pi], the authors show that 
on T 2 the global minimizer of the nonlocal isoperimetric problem (E^^-y with a = 0), agrees 
with the global minimizer of the isoperimetric problem, i.e., is given by ul, provided 7 > 0 is 
small. That is, the perimeter term dominates and the effect of the nonlocal perturbation via 
7 > 0 does not “kick-in” immediately whereas the above example shows that this is not the 
case for a > 0 . 


5. Concluding Remarks 


As noted in the introduction and as the example in Section [I] shows perhaps the most 
important feature of minimizing the energy E^ is that compared to the isoperimetric prob¬ 
lem the geometry of minimizing patterns can change significantly. Via its connection to the 
energy (Section §■ this reflects well the physical applications of adding nanoparti¬ 

cles into copolymer blends to change the morphology of pattern formation. Indeed, since the 


consideration of the energy ( 1 . 1 ) is to our knowledge the first mathematically rigorous study 


of nanoparticle/copolymer blends, this work also generates several directions for subjects of 
future studies. We will conclude by remarking on these directions. 

( 1 ) As mentioned before, depending on the choice of the penalizing measure p, the second 
term in E^ >(r can act as an attractive or repulsive term. For example, in two dimensions 
and small mass regime, by choosing the measure p distributed on disjoint small disks one 
can force the minimizer of the energy E A(i(T to “oscillate” rather than forming a larger disk 
which would be preferable in terms of minimizing the perimeter term (see Figure [5]). 

(2) Although the choice of the measure p provides a substantial freedom in forcing the mini- 
mizers of E^ jfr to form desired patterns, the minimizing patterns still exhibit some rigidity. 
In particular, regularity properties (Proposition [3TTj) and the criticality and stability condi¬ 
tions (Proposition |3.3| limit this freedom (see also the example in Section [I]). Adding the 
long-range interaction term between the phases controlled by 7 as in Remark 2.2 would 
enrich the possibilities for minimizing patterns. This would also provide further mathe¬ 


( 3 ) 


matical challenges in understanding the energy landscape of ( 2 . 8 ). 

Here we chose to fix the location of nanoparticles, hence, their distribution given by the 
measure p as the number of particles goes to infinity whereas their size approach zero. An 








SHARP INTERFACE OF NANOPARTICLE-POLYMER BLEND MODEL 


25 



Figure 5. A measure /i consisting of small disjoint supports (light gray blobs) 
might show a repulsive effect when a minimizer (striped regions) try to cover 
most of the support of p to reduce cost. 

interesting problem would be to analyze local and global minimizers of the energy E^ i0 - not 
only with respect to the phases, i.e., over u E BV( T n ; {±1}) with a fixed mass constraint 
m and fixed measure p, but also over the measures p E SP ac (T n ). 
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